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Abstract 

A solution to the classical field equations in the massless (l+l)-dimensional 
0(3) sigma model is found, which describes a multi-particle instanton-like tran- 
sition at high energy. In the limit of small number of initial particles, the 
number of final particles is shown to be also small, and the probability of the 
transition is suppressed by exp(— 2Sq), where Sq is the instanton action. This 
solution, however, does not correspond to the maximum transition probability 
among all states with given number of incoming particles and energy. Unless 
the limit g 2 ni n iti a i — ► is exponentially sensitive to the structure of the initial 
state, our results imply that well above the sphaleron energy, the instanton- 
induced cross section becomes again suppressed by the instanton exponent, and 
the number of final paricles is again small. 



1 Introduction 



Despite considerable effort, the problem of the rate of instanton-like transitions at high 
energies is still not well understood even in weakly coupled theories. The quantitative 
study of this problem was initiated by the observation [JXJ] (see also ref.@) that in 
the leading order of the pertubation theory around the instanton, the partial cross 
sections rapidly grow with energy, while the total cross section is exponential in energy 
[Q, [|, ^ |J and naively violates the unitarity bound at the sphaleron energy scale. 

Corrections to the leading order can be divided into three groups. The first group 
of corrections is entirely due to final particles and is called soft-soft (the multiplicity in 
the final state is large, and the momentum of final particles is low). These corrections 
sum up into an exponent, exp[F(E / E ) / g 2 ], where E is the sphaleron energy, and 
the function F contains contributions of tree graphs only|| [7|, §]. Two other groups of 
corrections, namely hard-hard (that involve only initial hard particles) and hard-soft 
(that involve both initial and final particles) contain non-trivial contributions from 
loops ||, 10] and are harder to deal with. There are, however, indications [10, 11] that 
they also sum up to an exponent of the same form as the soft-soft corrections, so 
one expects that the total cross section has the following form in the regime g —>■ 0, 
E/E = fixed, 

The value of F/g 2 at E = is of course (—2 So) where So is the instanton action; the 
behavior of F(E/Eq) at E <C Eq is calculable by the perturbation theory around the 
instanton in most models including the standard electroweak theory. 

To evaluate the total cross section at energies comparable to the sphaleron energy, 
one has to invent a new calculation technique. The exponential form of eq. QT|) suggests 
that there may exist semiclassical-type methods for calculating the total cross section 
at all energies. Clearly, such methods should involve a special treatment of the non- 
semiclassical initial state consisting of only two hard particles, as opposed to the 
multiparticle final state whose semiclassical treatment may be natural. 

One way to gain some insight into the problem of instanton transitions from 
two particle initial state is to consider the transitions from multiparticle states with 
n = v I g 2 initial particles, where v is small enough but fixed in the limit g 2 — > 
(and E / Eq also fixed) [0, [13|1 . In that case the problem of loop effects does not exist, 
and the transition probability from a given initial state can be evaluated, at least in 



principle, by entirely semiclassical technique outlined in refs.[12, 13, O] (see Sect. 2 
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This probability indeed has the exponential form, eq.([3]), where the function F now 
depends on the initial state and may be obtained by solving classical field equations 
with certain boundary conditions on a particular contour in complex time plane. An 
important issue in this procedure is the structure of singularities of the corresponding 
solution in the complex time plane. 

It may be quite natural that the exponent for the probability of the transition 
from states with vjg 2 particles is smooth in the limit v — > (after the limit g 2 — ► 
is taken) and that this limit is independent of the precise choice of the initial states 
(e.g., of the distribution of energy E over vjg 2 particles). In that case the limiting 
value F(E/E , v — > 0) would provide an upper bound to, and may even coincide with, 
the exponent for the two particle cross section. 

This assumption of the independence from the limiting value F[y — > 0) on the 
precise structure of initial states with vjg 2 particles, though plausible, is by no means 
justified. It may be more reasonable, therefore, to consider the "microcanonical" 
probability [T^, [T3, 13 1, i.e., the probability of the transition from a microcanonical 
ensemble at given center of mass energy E and the number of particles n = vjg 2 . 
Equivalently, one may study maximum transition probability among all initial states 
of given energy and number of particles. This probability, which is also calculable by 
semiclassical methods at finite v, almost by definition gives upper bound for the two 
particle cross section and is even more likely to coincide with the latter in the limit 
i/-»0. 

Substantial support to this line of reasoning has been recently provided by pertur- 
bative calculations about the instanton to high orders ((E / E ) 16 ^ 3 in the electroweak 
theory) RI5] . The conclusion of ref.fJTJJ is that in the perturbation theory about a 
single instanton, the process of transition of n\ = V\j ' g 2 left-moving and ri2 = 
right-moving W-bosons into n W-bosons goes smoothly to the 2 — > n process as 

Vt,V2 -> 0. 

To solve the relevant classical boundary value problem at generic values of E/Eq 
and v = ng 2 , one almost certainly should make use of numerical methods. However, 
it is desirable to find some analytic solutions that would provide benchmarks for 
numerical calculations and also indicate some properties of the instanton-like cross 
sections in limiting cases. At low energies, the approximation analogous to the dilute 
instanton gas is valid, and the corresponding solution for 2d Abelian Higgs model is 
given in ref.jli]]. The main purpose of this paper is to present a solution in the 2d 
massless 0(3) sigma model that describes an instanton-like transition in this model. 
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One may view this model as an extreme high energy limit, E/Eq 3> 1, of a massive 
sigma model |T6|, |i~7t jlSj , the latter sharing most of the relevant properties of the 
standard electroweak theory. 

A peculiar feature of the massless sigma model that made it possible to find 
an interesting solution is its conformal symmetry However, precisely this property 
prevents us to make contact to the low energy limit, E/Eq <C 1, of the massive 
theory where the transition rate is known from the perturbation theory about the 
instanton [17], |18 |. Furthermore, although our solution describes some instanton-like 
process, it does not correspond to any maximum (or "microcanonical") probability 
in the sense of refs.0, |13| . So, the relevance of our solution to the two particle 
instanton cross section relies on the assumption of the independence from the limit 
of small number of initial particles on the precise structure of the initial state. 

This paper is organized as follows. In Sect. 2 we review the coherent state for- 
malism in application to the study of multiparticle transitions at fixed energy. We 
derive the corresponding boundary value problem for the evaluation of the exponent 
for the total probability of a transition from a given initial state to all possible final 
states. We show that the sum over final states is in fact saturated by a single coherent 
state whose properties are determined by the classical solution. Inversely, for a given 
solution to the field equations with appropriate singularity structure in complex time 
plane, one can find the initial (and also final) state of the transition this solution 
corresponds to. We also point out a specific property of solutions that is required for 
their interpretation in terms of some "microcanonical" probability. 

Before turning to the actual solution, we discuss in Sect. 3 the peculiarity of asymp- 
totic states in (1+1) dimensional massless sigma model. Namely, asymptotic fields 
moving along the light "cone" x = t and x = — t are, in general, non-linear in terms of 
the original variables. We describe the change of variables that linearizes the asymp- 
totic fields (at least in the semiclassical approximation), thus making possible their 
interpretation in terms of particles. 

In Sect. 4 we obtain the solution explicitly and show that it has the desired singu- 
larity structure. 

Sect. 5 is devoted to the interpretation of the solution in terms of an instanton- 
like process. Specifically, we consider the case of small number of initial particles, 
v = g 2 n <C 1. We find that the number of final particles is also small, g 2 rifi na i <C 1 
(but non-analytic in u), and g 2 nfi na i vanishes in the limit v — > 0. The transition 
probability in this limit tends to the instanton value. Under the assumption of the 
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independence from the limit of small number of initial particles on the structure of 
initial state, this implies that the scattering at very high energies again becomes 
suppressed precisely by the instanton exponent. We show, however, that our solution 
cannot be interpreted as the configuration corresponding to the maximal transition 
probability among initial states with given number of initial particles; this makes the 
indication on the suppression of the instanton-like cross section at extremely high 
energies by precisely the instanton exponent less clear. 
Sect. 6 contains concluding remarks. 

2 Classical solutions and multiparticle scattering 
at fixed energy 



In the spirit of refs. [ 12 , [L3| , let us consider a sum 



*(W}) = E I «M I sp e I K» I 2 , (2) 

b 

where Pe is a projector onto a subspace of definite energy E, S is the instanton 
5-matrix and | {ak}), | {^k}) are coherent states. This sum is the probability 
of a transition at fixed energy E from some coherent state | {ctk}) (projected onto 
this energy) to all possible final states. We consider the system in its center of mass 
frame and will not discuss in what follows the constraints due to the total spatial 
momentum. 

Making use of the functional integral representation for the S-matrix elements, 
one writes the transition probability in the following integral form[^, |T2| 



dbldb k d^'d(j){x)d(l)'(x) exp ( - J dkb^ - iE(£ - £') + 

Bifae^tfa) + B f (&,<f> f ) +iS{<t>) + fl i *(a£e- iw **',#) + B} fa, fy) - iS (</>')), (3) 

where <ft is the condensed notation for the (real) bosonic fields of the theory, S((f)) is 
the action, B^ and Bf stand for boundary terms and are merely logarithms of the 
wave functions of coherent states in the representation, 

B l (a k ,(p l ) = -- j dk a k a_ k e- 2 ^ fcTl - -J dfcw k i (fc)0 i (-fe) 
+ J dky/2^e~^ kTi a^i{k) 



5 



Bf(bl 0/) = -\ Jdk btb\e 2 ^ T f - l -J dkukMQM-Q 

+ JdkV2^e luJkT fbl(j) f (-k) . (4) 

In the above formulae, Tj and Tf are initial and final times, fa, 0/ = 0(Tj), <p(Tf) and 
the limit T — > — oo, Tf — > +00 is assumed. 

At E ~ 1 / <7 2 and ~ 1/g, that corresponds to the regime 

g 2 ^ , v = ng 2 = fixed 

discussed throughout this paper, the integral in eq. (||) can be evaluated in the saddle- 
point approximation. Varying the exponent in eq.(|3]) with respect to 0, 0' at T < 
t < Tf, we obtain that the saddle point configuration for 0, 0' should be a solution 
to the field equations 5S/5(p = 5S/5(p' = 0, while variation with respect to 6£, fek 
and boundary values of and 0' leads to six equations which relate the boundary 
values of 0(x), 0'(x) and their time derivatives to a^, a£, 6 k , &£. In particular, the 
integration over b, b* , 0/ and faf results in the boundary conditions 

<j>'(x) = 00), fa(x) = 0(x) (5) 

at t — > +00, and 

%(Jfe) - cj k 0/(fc) + V^\b\e iuJkT f = 

- %(fe) - uj k (f> f {k) + v^&ke"^^ = (6) 

Eq.(|5|) implies that 

0'(x) = (f)(x) 

everywhere in space-time, i.e., we are dealing with a single solution to the field equa- 
tions. Eqs.(^) determine the saddle point values of the coherent state variables cor- 
responding to the final state. 

The saddle point equations at the initial time are 

ifa(k) - u^fa(k) + V^a\e M - luJke = 

- ifa(k) - u w fa(k) + V2^a k e~ M+i ^ = 0, (7) 

Notice that we keep the initial state, i.e., the amplitudes and a£, fixed, so that 
eqs. (0) are to be considered as the boundary conditions for the field at large negative 
time. 
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Finally, the variation with respect to £ and £' ( in fact, the integral over (£ + £') is 
trivial due to the time translation symmetry ) relates the energy E to other parameters 
of the solution. This relation is very simple and will be explicit later. 

The above boundary conditions can be converted into more convenient form. First, 
we assume that at the initial and final times T, and Ty, the field <fi becomes free and 
can be written as a superposition of plane waves. In particular, at large positive time 

0(M) = ^^(Ske-* 1 ** + g-^), t^+oo. 

Substituting this asymptotics into eq.@ we obtain 

h = Sht, bl = g k . (8) 

Notice that at this point we cannot insist that the saddle point field is real, so in 
general g k need not be complex conjugate to ^it- 
Second, the dependence of the integrand in eq.(^) on (£ + £') is trivial, so we can 
set £' = — £. Furthermore, the real part of the saddle point value of £ can be removed 
by the time translation, and we can choose 

i = iT 

Consider now the contour in complex time plane shown in fig. 1. On this contour, the 
field at early time, 

t = iT + t', t' -> -oo 

has the following form, 

<f>(k, t) = -^(/ k e-^*' + /_ k e^') (9) 

where / k need not be complex conjugate to / k . The second of eqs.(|7|) gives then 

/k = a k (10) 

Analogously, the first boundary condition in eq.(^) is conveniently formulated at 
Im t = —T, where the field is 

0(M") = -^={Ke-^' + /L k e*"* t ") 

t = -iT + t", t" -> -oo 
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The boundary condition reads 

h k = al 

Since a k ans a k are complex conjugate to each other, the latter condition is equivalent 
to the requirement that the solution <p{x) is real on the real (minkowskian) part of 
the contour, C in fig. 1, provided the solution is unique. This means, in particular, 
that the saddle point values of 6 k and 5 k are complex conjugate to each other, i.e., 
the sum over final states in eq.(§) is saturated by a single coherent state. 

To summarize, the transition probability cx({a k }) is determined, in the semiclassi- 
cal approximation, by a solution to the classical field equations on the contour of fig.l, 
obeying the following boundary conditions: (i) the field must be real in minkowskian 
part of the contour; (ii) the negative frequency (" anti-Feynman" ) part of the field at 
large negative time (part A in fig.l) must be determined by the amplitudes a k via 
eqs.(P),(|ToD. The final state is then given by the asymptotics of the solution at large 
positive time by eq.(||). 

Finally, one can see that the saddle point value of T is determined by the require- 
ment that the classical energy of the solution is equal to E. 

A few remarks are in order. First, the solution need not be real on the euclidean 
part and/or on the part (AB) of the contour of fig.l. Furthermore, the solutions 
real everywhere on this contour do not correspond to the case when the number of 
initial particles is small |19|, |T2| . So, most interesting solutions even in theories with 
hermitean fields are complex. Second, interesting solutions, though smooth at our 
contour, should have singularities in the complex time plane between the real axis 
and the line {AB) , otherwise they would correspond to classically allowed scattering. 

Finally, instead of considering transitions from a given initial state, on may study 
the "microcanonical" probability [|T^ , (0|, 

a E , D = J2\(f\SP E PD\t)\ 2 (11) 

where the sum runs over all final and initial states and Pjj is the projector onto an 
eigenspace of some operator 

D = J dkA k AlA k 

with real A k (an obvious choice would be A k = 1 that would correspond to fixing the 
number of initial particles). Here A^ and Ak are creation and annihilation operators 
for initial particles. It follows from refs.[12|, [13], [14] that the sum in eq.flTTD is also 
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saturated by a single initial state which therefore has the maximum transition prob- 
ability among the states with given energy E and the value of D. Instead of eq . (|10|) , 
the boundary condition for the corresponding saddle point field would be 

/k = e" eAfc / k (12) 

in the region (A) of the contour of fig.l, where 9 is some real constant and f£ is 
complex conjugate to / k . Notice that in this case the ratio of positive ("Feynman") 
and negative (" ant i- Feynman" ) frequency parts of the saddle point field at t = iT + t', 
t' — > — oo, should be real and positive. 

Let us now come back to the discussion of the transitions from a fixed initial 
state. The above arguments show that every solution to the field equations, that 
is real on real time axis and has an appropriate singularity structure in complex 
time plane, corresponds to such a transition from some initial state. The amplitudes 
characterising the initial state, {a k }, are determined by the negative frequency part 
of the solution in the region A of the contour via eq.fllOD, while the most probable 
final state is determined by the behavior of the solution at large real time by eq.(§). 
There is, however, an apparent ambiguity in the determination of the initial state of 
the process described by a given solution. Namely, one can shift the line AB of this 
contour upwards or downwards (until this line reaches the singularity) thus changing 
the asymptotics of the field along the contour at negative Re t. This would induce 
a change in the amplitudes a k of the initial state. In other words, the choice of the 
value of T for a given solution has apparent arbitrariness, that translates into the 
apparent arbitrariness of the initial state. 

Let us see that this ambiguity is in fact irrelevant. The key observation is that 
the projections of the coherent states of the form | a^^) with various rj onto the 
subspace of the hamiltonian H with energy E differ only by normalization: 

Pe I a k e^> = / dte l(H ~ E)t | a k e^) = 

= J dte~ iEt | a k e iLUkt+ " k ' n ) = J dt 'e~ iE{t >+ir >) | a k e^'') = e Eri P E | a k ) 

Suppose now that there are no singularities of the field in the region T < Im t < T'. 
In this case the negative frequency components / k and of the solution on the lines 
Im t = T and Im t = T', respectively, are related to each other as follows, 

fL = he^\ r) = T'-T 
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The same relation holds for the amplitudes characterising initial states, so their pro- 
jections onto the subspace of fixed energy E are essentially the same. 

To find the most convenient choice of the contour (i.e., the most convenient value 
of the parameter T), let us first fix T in an arbitrary way and calculate the average 
of some bilinear operator 

= J dkF(k)AlA k 

over the state Pe | cik) determined by the solution on the contour that starts at 
Im t = T. We have 

(a k | P E OP E | a k ) = 
J dtdt'e- iE{t ~ t '\a k e luJkt ' \ J dkF{k)A{A k \ a k e iWkt ) = 

J dtdt'exp (-iE(t -t') + J dkf^e^^'A J dkF^flf^^. (13) 

where is complex conjugate to f k (notice that f£ are not equal to the positive 
frequency components of the solution, f^). The integral can be evaluated in the 
saddle-point approximation. The result is 

(O) = J dkF{k)fU^ T \ 

where T is detrmined by the following relation, 

E = J dku k fZf k e-^ To , (14) 

(the factor exp{ET + f dkflfae-^ ) that appears in eq.([L3|) at the saddle-point is 
merely the norm of the state Pe \ Ok)). If instead of a k = f k we used a k = f k e~ u)kT °/ 2 
(so that Pe \ ctk) and Pe \ «k) are the same states with different normalizations), 
then the average of O would have a simple form, 



(O) = J dkF(k)ala w . 



Note that the average energy of the coherent state | «k) is equal to E (eq.(|nj)). If 
there are no singularities in the region T — To/2 < Im t < T, the positive frequency 
components of the solution calculated at t = i(T — To/2) + t', t' — > — oo are exactly 
equal to «k- So, in some cases it is possible to choose T in the most convenient way, 
namely, to satisfy the condition 



E = J rfk^k/k/k- 
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With this choice of the contour, the averages of operators over the initial state (say, 
the number of particles) are expressed through (ak, a k ) = (/k, /£) in a standard way. 

So, every solution that is real on the real time axis and has correct singularity 
structure in the complex time plane (so that one can really draw a contour of fig.l 
at which the solution is non-singular) can be considered as a configuration that dom- 
inates some process at fixed energy The semiclassical probability of this process 
can be found by substituing the saddle point values of the integration variables into 
eq.(|3|) and dividing the result by the norm of the state Pg | a k ), which is equal to 
exp(J dka^a k ) for the contour with "correct" T. One obtains for the semiclassical 
probability 

a = exp (W - 21m S(<f>) + J dk(/ k - /£)/ k ) . (15) 

If, furthermore, at large negative time (part A of the contour) positive and negative 
frequency components for a given solution are complex conjugate to each other up to 
a real and positive factor for all k, then this solution can be interpreted as describing 
some microcanonical probability: indeed, in that case one can choose a suitable form 
of Ak to satisfy eq.([E|). We shall see that this property is absent for the sigma model 
solution of Sect. 4. 

To apply this formalism to the (1+1) dimensional sigma model we have, however, 
to be able to describe its asymptotic states. 



3 Asymptotic states in 0(3) sigma model 

Consider the massless non-linear 0(3) sigma model with the action 

1 r » ,o 2 r , d,,wd,,w 



s = —l d^nf = - 2 / f 

2g z J g z J (1 + ww) z 

n 2 = l 

where w and w are fields in the projective representation, 

w + w 1 w — w 1 — WW 

ni = : , _ , n 2 = -— — , n 3 = — — . 16 

1 + WW % 1 + WW 1 + WW 

A peculiar feature of this model is that the fields cannot be considered free in asymp- 
totic regions t — > ±oo. Indeed, any field configuration of the form w = w(x + ) or 
w = w(x-) where x + = x + t, x_ = x — t, satisfies the field equations. Its amplitude 



11 



need not be small, and the non-linearity persists even at large positive and/or nega- 
tive t. The interpretation of classical solutions in terms of particles requires a change 
of variables. 

Let us consider only such solutions that in the limit t — > oo split into two non- 
linear wave packets moving along two parts of the light "cone" x — t and x = —t, 

w(x, t) = w + (x + ) + W-(X-). 

At large time the two wave packets do not interact, so it is natural that left moving 
and right moving waves can be considered separately. Let, for example, w = w(x ) 
be right moving quantum field. Let us examine the generator of the translation along 
the X- direction, 

P.-Xftx- -* 8 -" 



g 2 J (1 + to) 2 ' 
It has the following property, 

[P-,w(x-)] = id^w(x_), 

[P-,w(x-)] — id-w(x-), 

where we have used the notation <9_ = d/dx-. Let us now define new fields <&(#_) and 
by the following relations, which are unambiguous at least in the semiclassical 

limit, 

g 1 + ww g 1 + ww 

Then the generator P_ becomes 

p_ = 2 J d 2 xd^d^, 

while its commutators with $ and $ have the same functional form as with w and w, 

[P_,$(x-)] =id-<S>(x^), 

[P_,|.(x_)] =id_$(x_). 

To ensure these commutation ralations, it is sufficient to require the following com- 
mutation ralations for $ and $, 

[$(x-)My-)] = i<x--v-)- (is) 
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In the Fourier representation, 



/■OO rlh 

*(*-)=/ -7-7=^-+^-), 



Jo V27TV2A; 



the commutation relations (0) imply that the operators a k , 6r k , /3 k , $ k can be inter- 
preted as annihilation and creation operators for particles and anti-particles, respec- 
tively, 

[a k ,aU=6(k-k'), k Jl} = 6(k-k'), 

and other commutators vanish. Thus, the change of variables, eq. (|TTD , enables us to 
interprete the outgoing right moving fields in terms of particles and antiparticles of 
spatial momentum k > 0. Quantizing the left-moving waves we describe particles and 
antiparticles with negative momenta in the same way. Note that our interpretation 
is essentially based on the separation of two wave packets in the asymptotic region 
(eq-©)- 

Once we have defined annihilation and creation operators, we can construct co- 
herent states in the standard way, 

a k | a k ) = a k \ a k ), (3 k \ (3 k ) = (3 k \ fi k ), 

and make use of the coherent-state formalism for calculating transition probabilities 
at fixed energy in the sigma model. 



4 The classical solution 

Our approach to obtain the classical solutions that are real (in the sense that w is 
complex conjugate to w) in Minkowski time and have singularities in the complex 
time plane is to start from euclidean space-time (x, r) and search for real solutions 
that have a turning point, d T w(x) = for all x at r = 0. The reality of the solutions 
in the euclidean domain and the existence of the turning point then guarantee that 
their analytical continuation to the complex time plane is real at Minkowski time. 
Since we are interested only in a finite region of euclidean time, < t < T, the 
solutions may (and, in fact, should) have singularities outside this region. 
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Since the action of the sigma model is conformally invariant, the field equations 
do not change under the following conformal mapping, 



z — i x + ir — i 



X 



IT 



(19) 



Z + % X + IT + I 

that tranforms the upper half-plane into a disc of unit radius, centered at z' = 0, in 
which we solve the field equations. [There is no intrinsic mass scale in the model at 
the classical level, so we do not introduce a length scale into the solution. We will 
be interested in dimensionless quantities only.] The existence of the turning point in 
the original variables is equivalent, in terms of the new variables, to the boundary 
condition that the normal derivative of w at the circle of unit radius that bounds the 
disc, is equal to zero. If we introduce the polar coordinates p and (p on the z'-plane, 



P 



v + (t-i; 
x 2 + (t + r 



1/2 




1/2 



X + lY + T l 



then the boundary condition becomes d p w \ p =i- 
We search for 0(2) symmetric solutions, 



0. 



w 



■up 



The action for these configurations has the form 

S 



4?r f 1 , 1 

f ' dp TTTW 




+ 



P 



Making the change of variablesfQ / = tan(-^/4), p = exp(£) one obtains 



7r r°° 



1 fdf 

2 [dC, 



+ (1 — cosip) 



The problem is now reduced to the mechanical problem of a particle moving in the 
sin-Gordon potential U(Q = —(1 — cosip) with the initial condition dip/d( = at 
C = (Fig. 2). There exists a set of solutions depending on a single parameter ip 
which is the initial coordinate of the particle. It is obvious that the solution ip(() is 
periodic. The oscillatory behavior of V'(C) means that the point p = (£ = — oo) is an 
essential singularity. In the initial z-plane this essential singularity appears at x = 0, 
r = 1. The existence of essential singularities in the euclidean time axis may, in fact, 
have been expected for configurations describing processes in which the number of 
initial and final particles are not equal to each other [|i~4|j . 
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The explicit form of the solution is 

ij) — 2 arccosf— k sn(( — K, k )], (20) 
where k = cos(V>o/2), sn is the Jacobi elliptic sine, and 

K -- 



f/ 2 dx 



\Jl-kl sin 2 x 

is a quarter of the period of function ip(C)- The fields w and -u) are then, 



w 



l + k sn((-K) \ 1/2 _ f l + k sn((-K) \ 1/2 



where £ and 99 must be understood as functions of x and r, 

1 x 2 + (r + l) 2 i (x + Q 2 + t 2 

C = lnp= 2 ln x 2 + (T-ir ^ = 2 ln (x-^ + r 2 - (22) 

Now we can continue analytically the solution to complex time plane. This is 
done by setting r = —it where t will finally run along the contour of fig. 1. Eq.(p2[) 
transforms into 

j., x 1, x 2 + (l-itf i, (x + i) 2 -t 2 

2 a;" 1 + (1 + zt)^ 2 (x — «) z — t z 

When t is complex, f and C are, in general, also complex, so the fields w and w are 
not complex conjugate to each other. 

Let us now turn to the discussion of the singularities of our solution in complex 
time. There exist essential singularities at 

t = ±x + i 

(we concentrate on the region Im t > 0). Their positions make a "cross" in the three 
dimensional space (Re t, Im t,x), as shown in fig. 3. Clearly, we should choose T < 1 
for the contour of fig.l not to intersect this cross. 

One can see that apart from these essential singularities, the fields w and w have 
poles at complex t. However, at these points the original fields n and other dynamical 
variables (like energy density) are regular. Therefore, these poles have no physical 
significance. 

On the other hand, the fields n and other dynamical quantities are singular at the 
points where 

1 + ww = 
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i.e., where the transformation (plj) is singular. At these points the dynamical variables 
have poles, and the positions of these poles, t = t(x), make (non-planar) curves in 
the space (Re t, Im t, ; clS shown schematically in fig. 3. 
We will be interested mostly in the case 

< 1 

which, as will be shown in Sect. 5, corresponds to small number of initial particles. 
Let us describe the positions of the physical singularities in that case. Making use of 
the explicit form of the solution, we write 

9 

1 + WW 



l-k sn{(-K) ' 

which is zero when sn(( — K) is infinite, or ( = (2n + 1)K + (2m + l)iK', where 

r' 2 dx 

K 



1 — k 2 sin 2 x 
k\ = \ — k\ = sin 2 

At small ^o, the equation ((x,t) = (2n + 1)K + (2m + l)iK', where ((x,t) is given 
by eq.(^), has solutions only at m = — 1 and n > 0. The singularity curves 

t = t n (x), 72 = 0,1,2,... 

have Im t close to 1, and the larger is n, the closer is the curve t = t n (x) to the cross 
t = ±x + i (two of these curves are shown in fig.3; our solution is symmetric under the 
transformation t —>■ —t, so there are also singularity curves at Re t > 0, not shown in 
fig.3). Asymptotically, these curves run parallel to the light "cone", and the behavior 
of t n (x) at large | x \ is 




t n = — \ x \ +i \ 1 — 2 — + • ■ • \ x — > oo. 



The integration contour of fig.l can be imagined as a ladder step of height T in 
the three-dimensional picture of fig.3. By varying T we can place the upper part of 
the step, (Re t < 0, Im t = T), above different number of the singularity curves. To 
avoid crossing the singularities, one should in general consider curved surfaces, not 
just ladder steps, but this is unnecessary in our case for the surfaces with exactly one 
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singularity line under them. In what follows we will consider only such surfaces, i.e., 
we will choose 

1 - k 2 j8 < T(oo) < 1 - kf/2 11 . 

This choice corresponds to one-instanton transitions, while surfaces with more than 
one singularity curves under them would most likely correspond to multi-instanton 
processes. 

The solution is explicitly real (in the sense that w is complex conjugate to w) in 
Minkowski time. This is clear from the following equivalent form of the solution, 

and the observation that at real t, the variables ( and tp determined by eq. (^3|) are 
pure imaginary and real, respectively. 

To calculate the imaginary part of the action and the topological charge Q, we 
notice that for every x, the integral 

b [X) ~ Jc (1 + w(x, t)w(x, t)f ' [ b) 

on the integration contour of fig.l is a sum of the same integral over the real time axis, 
and the contribution from the pole 1 + ww = 0. On the real time axis the integral 
is real, so the imaginary part comes from the pole. [The integral for the topological 
charge over the real time axis is zero because of the invariance of the solution under 
t — > —t.} Then the calculation is straightforward, and the result is 

47T 

? 



f 47T 

Im S = / dxSix) = — 

J Q 1 

In the same way we find 



Q = l 

Thus, the action is exactly equal to the instanton action for all ipo, and our solution 
indeed corresponds to a one-instanton transition. 
At large | x \ and fixed t, the fields tend to 

w = w = w Q 

where 

w = tan(W4) (26) 
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Therefore, our configuration is a fluctuation about a vacuum with non-vanishing 
expectation values of w and w, i.e., about a vacuum with (0, 0, 1). It is convenient 
to perform an 0(3) rotation in the n-space so that the vacuum expectation value 
becomes n = (0, 0, 1). In terms of the projective variables this rotation reads 

w = , w = . (27) 

1 + w w 1 + w w 

The fields w' amd w' are then the fluctuations about their zero vacuum expectation 
values. 

Let us discuss the solution in the case ipo <C 1 in more detail. In that case the 
singularity line is near the plane Im t = 1. Well below this plane, we make use of the 
following approximation for dn(i(,ki) at small k\ |p^ |, 



dn«,A; 1 ) = l + ^sinh 2 C, (2? 



so that to the leading order in ip we have 



, 1 1 

W = -IWq ; + 



x + t + i x — t + % 

w' = tw ( J. . + , (29) 

\x + 1 — % x — t — ij 

Wo = ^ipo 

In the euclidean domain, t = ir, this field configuration is that of widely separated 
instanton and anti-instanton of sizes Wq <C 1, sitting at x = and r = +1 and 
t = — 1, respectively. This property may have been expected on the basis of the 
results of ref. |14]]. At (1 — r) ~ Wq we cannot, however, use the approximation (p8[) 
and the solution differs from the instanton - anti-instanton pair considerably. 

To end up this section, let us point out that our solution is very similar to the 



Luscher-Schechter solution in (3 + 1) pure Yang-Mills theory|21, 22]. In fact, we could 
generate more solutions in the sigma model by allowing the integration constant in 
eq.(|20|) to be (— K + i( ) instead of (— K). For any real ( these configurations 



are real in Minkowski time. In analogy to the Luscher-Schechter solutions, they 



possess a peculiar feature, discovered in ref.^3j for (3 + l)d Yang-Mills theory, that 
their topological number, evaluated in Minkowski time, is non-integer in the sense of 
ref.p3[- We will not discuss this property in the present paper and concentrate on 
the case Co = 0. 
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5 The semiclassical process 



To interprete our solution in terms of a semiclassical process, we have to evaluate the 
asymptotics of the fields $ and $ which are constructed according to eq.([L7]) from 
w' and w' defined by eq.(^7|). As discussed in Sect. 2, their Fourier components at 
t — > +00 are the amplitudes of the final coherent state, while their anti-Feynman 
parts at t = iT + 1', t' — > —00 determine the amplitudes of the initial coherent state. 
We consider the case ipo "C 1 that corresponds to the small number of initial particles. 



5.1 Final state 

We begin with the discussion of the final state. To the leading order in ip (or in w , 
see eq. (|26|) ), the field at real time is given by eq. (|29|) . At large t it is a sum of waves 
moving left and right along the light "cone". This field is weak, and to the leading 
order in w 

$ = V2w'/g, <5 = V2w'/g 

The Fourier components are easy to calculate, 

ttfc = at = — ^ woVke~ k , (3k = (31 = 0, for k > 0, 
9 

/ 87T 



a k = at = 0, (3 k = (31 = w V-ke k , for k < 0. (30) 

9 

We see that the final state contains particles moving right and anti-particles moving 
left. The energy and total number of particles are 

/oo 47T 
dk\k\ (a* k a k + (3l(3 k ) = —w 2 , (31) 
-00 Q 

/oo 47T 
dk(a* k a k + (3* k (3 k ) = -^w\. (32) 
-00 g 

The typical momentum of a final particle is 

k final ~ 1 (33) 

Of course, the latter estimate has no direct physical significance, but it will be useful 
for the comparison to the typical momentum of initial particles. 
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5.2 Initial state 

The calculation of the Fourier components of the fields $ and $ for the initial state 
is less straightforward. At w -C 1, the fields w and w are large for Im t close to 1. 
So, we really have to use the formalism developed in Sect. 3. 

Let us consider the regime t = iT + t', if — > — oo and concentrate on the left part 
of the light cone where x is close to t' . In this region 

1 x — % 

((x,t) =i<p= -In — , 

2 x_ + 1 

where 

X- = x — i! — iT 

At X- far from the singularity, we may still use the approximation fl28|) and obtain 

w = w X ~ . , w = w , , (34) 

Notice that at T 7^ the fields w and w are not comlex-conjugate to each other. 

However, the approximation (p8|) is not valid when ( is close to the singularity 
(K — iK') which happens if (1 — T) is of order or less and x_ near —i. Let us 
evaluate w and w for these values of (. Since £ is close to (K — iK r ), the most suitable 
(and equivalent) form of the solution is as follows, 



.„ l + m(C-K + iK',ko) V /2 
-1 + sn(C - if + zlf', V 



where we made use of the transfomation formulae for elliptic functions (see ref . [[2~4|j ) . 
At small tpo, the leading-order approximation for the elliptic sine is sn u = tanhu, 
and we obtain 

-iw _ wo 

w = — , w = — , 35 

x_ + 1 2 



This is consistent with eq.fl34) when x_ is close to —i, which means that the "naive" 
formulae also work in the region near the singularity. However, to calculate $ and $ 
we have to find the derivatives of w and w with respect to a;_. For d_w we can use 
either eq.([34D or eq.(^), but the leading-order expression obtained from eq.(^) for 
d-W vanishes, so we have to make use of the second-order approximation! 



sn(w, ko) = 1 + — ) tanh 
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from which we obtain for X- near — z, 

w o ( \ w o i x - + i i w o 1 i , ■^^ 

w = — 1 + — In 1 h -(ac_ + 1) , 

2 V 2 -2z 8 x_ +i 2 V 7 ' 

and 

_ _ z'^o + * - ^o/ 2 ) 2 
cLw = — t — . 

4 (x_+z) 2 

These expressions do not coincide with eq. ([34]) , so we have two approximate formulae 
for d-W in two different regions, x_ + z ~ w 2 . and | x_ + i |^> zu 2 . One can construct 
a function that approximates cLw uniformly in the two regions, namely, 

d w = -i ( x - + ^ ~ ^ w o/ 2 ) 2 
(x_ — i) 2 (x- + i) 2 

This expression can now be used for calculating $ and $. After performing the 
rotation, eq.(^7|), we get 

_ , \/2 <9_k/ \/2 1 + w w cLw 



g 1 + u>W g 1 + u>ow 1 + 

i^2 



^ (x_ + z — iw1/2){x- + z — zwq) 



2\ ' 



<9_$ 



re' 



j 1 + 

(x_ + i — zWq/2)(o;_ + z — zw 2 ) 
(x_ — z') 2 (x_ + z') 2 

The calculation of the Fourier components is now straightforward. We find for 
the anti-Feynman parts 

OLk = 0, 

P k = -\^w 5 o^7t ( k -X) ( 36 ) 



^V8 °v^V 

and for the Feynman parts 



w. 



a k = -w Vke (2 e)k , 

9 



fa = -7327—^ ( e -K 2 /2-,)/c _ e -H-e)t\ (37) 



where 

e = 1-T« 1 
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is yet an arbitrary parameter in the interval (wq/32, Wq/2). Eqs.(|36|),(|3~7D are valid 
for the right-moving waves (k > 0), the expressions for the left-moving waves (k < 0) 
are completely analogous. 

As pointed out in Sect. 2, the contour in the complex time plane, i.e, the value of 
the parameter e, should be chosen in such a way that the initial coherent state has 
the same average energy as the final one. In our case this means 

roo 47T 

dk\k\ (a* k a k + f3* k f3 k ) = —^w\, (38) 

g 

Recall that a* k and (3% are complex congugate to a k and /3 k , respectively, and are not 
equal to the Feynman parts of the field. Eq . (|38|) gives 

Eq. (|36|) then implies that the typical momentum of an incoming particle is 

; -8/3 
^initial ~ ^0 

which is much larger than the typical momentum of an outgoing particle, eq. fl3"3|) . 
The number of particles in the initial state is much smaller than the number of final 
particles, 

n = / dk(a* k a k + (3* k f3 k ) = -^w 1 ^ 3 . (39) 
J — oo g 

Eq. (|39| ) shows that the process with small number of incoming particles indeed cor- 
responds to the limit Wq — * 0. In this limit the number of outgoing particles is also 



small, see eq.(|32D, so that 

n fi na i = -^{g 2 n initial f 17 (40) 
T 

Notice that the number of final particles is non-analytic in v = g 2 ni n i t i a i- 

Finally, we note that our solution does not correspond to any "microcanonical 
process" in the sense of Sect. 2. Indeed, the ratio Pk/P k is not a positive definite 
function of momentum, which would be necessary for the interpretation of the solution 
in terms of the "microcanonical" probability. 



5.3 Transition probability 

The probability of the process described by our solution is calculated according to 
eq. (|T5|) . We have seen in Sect. 4 that the imaginary part of the action for our con- 
figuration is exactly equal to the instanton action. The parameter T at small u>o 
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IS 



T = 1 + OH 8/3 ) 



and eqs.(|3TD,(|3"6D,(|3~?D and give for various terms in the exponent 



2ET = — 



Jdk(a k a k + PM mtUal = (l-31n2 



and 

dk(a* k a k + (3* k (3 k ) initial = 0(w 1 q 13 ) 



So, we have for the transition probability 

/ Stt 127r(l-ln2) 2 14/3 \ 
a = CX P ( + p + O(w ) J 

Equivalently, the transition probability may be expressed through the number of the 
incoming particles, 

a = exp | + 12(1 - In 2) (^j ' {n imtial f 7 \ (41) 

Clearly, at small n initiah the transition probability is suppressed by the instanton 
exponent. 



6 Conclusions 

In this paper we have found an analytic solution to the classical field equations in the 
massless 0(3) sigma model in (1 + 1) dimensions. This solution describes a multi- 
particle instanton-like transition at high energy. Most properties of this solution agree 
with the earlier expectations concerning the structure of its singularities in the 
complex time plane and the behavior in the euclidean domain in the limit of small 
number of incoming particles. We have found, however, that this solution does not 
correspond to the maximum transition probability among all initial states with given 
number of paricles and given energy (or any "microcanonical" probability in general). 

Unless the limit v = g 2 n ini tiai — > is very sensitive to the details of the structure 
of the initial state (which seems rather unlikely), our results indicate that well above 
the sphaleron, the instanton-induced cross section again becomes suppressed by the 
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instanton exponent, and the number of outgoing particles is again much smaller than 
const/g 2 . However, on the basis of these results one may anticipate, for 2 — > any 
process, the non-analiticity in g 2 of both the number of outgoing particles and the 
exponent for the total cross section. Indeed, if we naively set n in iu a i ~ 1 in eqs.(f40"D 
and we would get 

1 

n final ~ 

and 

/ 8n const \ 

To obtain the solution analytically, we heavily made use of specific properties of 
the massless theory, in particular, its conformal invariance. Therefore, we were unable 
to make contact to the perturbative calculations about the instanton that are relevant 
at low energies. Hopefully, the complete solution of this problem may be possible by 
numerical methods. 

The authors are indebted to E.Farhi, V.Khoze, L.McLerran, A.Mueller, A.Ring- 
wald, M.Shaposhnikov, A.Vainshtein and M.Voloshin for discussions. P.T. thanks 
CERN for hospitality and the Eloisatron project for partial support, V.R. thanks 
TPI, University of Minnesota, for hospitality. 

References 

[1] A.Ringwald, Nucl.Phys. B330 (1990) 1. 
[2] O.Espinosa, Nucl.Phys. B334 (1990) 310. 

[3] L.McLerran, A.Vainshtein and M.B.Voloshin, Phys. Rev. D42 (1990) 171. 

[4] V.I.Zakharov, Classical Corrections to Instanton- Induced Amplitudes, preprint 
TPI-MINN-90/7-T, 1990; 
Nucl.Phys., B371 (1992) 637. 

[5] S.Yu.Khlebnikov, V.A.Rubakov and P.G.Tinyakov, Mod.Phys.Lett. A5 (1990) 
1983; 

Nucl.Phys. B350 (1991) 441. 
[6] M.Porrati, Nucl.Phys. B347 (1991) 371. 
[7] P.B.Arnold and M.P.Mattis, Phys. Rev. D42 (1990) 1738. 



24 



[8] L.Yaffe, Scattering Amplitudes in Instanton Background. In: M.Mattis and 
E.Mottola, eds., Proc. of the Santa Fe Workshop, World Scientific, Singapore, 
1990. 

[9] A.Mueller, Nucl.Phys. B348 (1991) 310. 
[10] A.Mueller, Nucl.Phys. B353 (1991) 44. 

[11] S.Yu.Khlebnikov and P.G.Tinyakov, Phys.Lett 269B (1991) 149; 
P.B.Arnold and M.P.Mattis, Phys.Rev. D44 (1991) 3650; 
M.P.Mattis, L.McLerran and L.G.Yaffe, Phys.Rev. D45 (1992) 4294. 

[12] V.A.Rubakov and P.G.Tinyakov, Phys. Lett. B279 (1992) 165. 

[13] P.G.Tinyakov, Phys. Lett. B284 (1992) 410. 

[14] V.A.Rubakov, D.T.Son and P.G.Tinyakov, Phys.Lett 287B (1992) 342. 

[15] A.Mueller, Comparing Two Particle and Multi-Particle Initiated Processes in the 
One Instanton Sector, preprint CU-TP-572, 1992. 

[16] E.Mottola and A.Wipf, Phys.Rev. D39 (1989) 588. 

[17] S.Yu.Khlebnikov, V.A.Rubakov and P.G.Tinyakov, Nucl.Phys. B347 (1990) 783. 

[18] N.Dorey, Multi- Instanton Valleys in the 0(3) Sigma-Model, preprint LA-UR-92- 
1380, 1992. 

[19] S.Yu.Khlebnikov, V.A.Rubakov and P.G.Tinyakov, Nucl.Phys. B367 (1991) 334. 
[20] J.Bossart and C.Weisendanger, Phys.Rev. D46 (1992) 1820. 
[21] M.Liischer, Phys.Lett 70B (1977) 321. 
[22] B.Schechter, Phys.Rev. D16 (1977) 3015. 

[23] E.Farhi, V.V.Khoze and R. Singleton, Jr., Minkowski Space Non-Abelian Classical 
Solutions with Non-Integer Topological Number, MIT preprint, 1992. 

[24] Handbook of Mathematical Functions, M.Abramowitz and I.Stegun, eds., Dover, 
N.-Y.,1964. 

I.S.Gradstein and I.M.Ryzhik, Table of Integrals, Series and Products, Academic 
Press, N.-Y., 1980. 



25 



B 



Im t 
T 



C 



Re t 



Fig.l 



26 




27 




28 



